All-electrical control of single ion spins in a semiconductor 
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We propose a method for all-electrical manipulation of single ion spins substituted into a semicon- 
ductor. Mn ions with a bound hole in GaAs form a natural example. Direct electrical manipulation 
of the ion spin is possible because electric fields manipulate the orbital wave function of the hole, and 
through the spin-orbit coupling the spin is reoriented as well. Coupling ion spins can be achieved 
using gates to control the size of the hole wave function. Coherent manipulation of ionic spins may 
find applications in high density storage and in scalable coherent or quantum information processing. 
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Observing magnetic resonance between different spin 
states of nuclei or electrons is a technique widely used 
in many imaging and spectroscopic applications. Sen- 
sitivity sufficient to measure the fluctuation of a sin- 
gle spin has been demonstrated using magnetic reso- 
nance force microscopy Q, noise spectroscopy 0, op- 
tical spectroscopy |3|, ]^ , scanning tunneling microscopy 
(STM) HH0, and quantum point contact conductiv- 
ity ||. Thepotential also exists to determine the spin 
state 0, 0, S • Controlling a single spin, in addition 
to monitoring it, is highly desirable for building future 
spin-based devices, is essential for quantum computa- 
tion, and should permit the direct exploration of fun- 
damental aspects of quantum dynamics in a solid state 
environment |l(J, 111! Il2l 1 1 3| . Proposed schemes to con- 
trol a single spin in a solid state environment rely ei- 
ther on magnetic resonance 0, or optical manipu- 
lation [lfj, [lj], , and there has been progress towards 
replacing ma gne tic control fields with the spin-orbit in- 
teraction [llll20l] or the exchange interaction 0, |53 • 

Ionic spin states in solids have several attractive char- 
acteristics for fundamental studies of spin dynamics and 
for spin-based devices. Every ion embedded in a solid is 
identical to every other such ion. Thus an ionic spin sys- 
tem can be as uniform as a nuclear spin system, but also 
can permit spin manipulation on short time scales as in a 
quantum dot spin system. Controlling ionic single spins 
without any magnetic fields, using techniques in which 
electric fields play the typical role of magnetic fields, may 
therefore provide a path to high-density scalable spin- 
based electronics. For example, the control of ionic spin 
states can be used to produce highly spin-selective and 
spin-dependent tunneling currents in nanoscale electrical 
devices, or to realize quantum computation. Manipula- 
tion of individual spins that are constituents of interact- 
ing spin clusters also opens up the capability to explore 
the fundamental dynamics of frustrated spin systems and 
other correlated spin systems. 

Here we propose an all-electrical scheme for ionic spin 
manipulation in which the role of magnetic fields in tradi- 
tional electron spin resonance (ESR) is replaced by elec- 
tric fields. In conventional ESR the energy splitting be- 



tween different spin states, and the couplings between 
them, are controlled by magnetic fields because an elec- 
tric field does not directly couple to the electron's spin. In 
a semiconductor crystal with tetrahedral symmetry and 
spin-orbit interaction (such as GaAs) a J = 1 ion spin 
(such as that of Mn in GaAs) will be triply degenerate, 
however the energy splittings and the couplings between 
these states depend linearly on the electric field strength, 
allowing rapid all-electrical control. Thus all operations 
performed with magnetic fields in traditional ESR, can 
be performed with electrical techniques. 

A specific proposed setup for manipulating a single ion 
spin is shown in Fig. ^ Tip-induced placement of single 
Mn ions substituted for Ga in a GaAs sample has been 
demonstrated experimentally [2^. Two gates are config- 
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FIG. 1: (color online) Proposed configuration for the electric 
resonances of a single Mn dopant in GaAs. A dc electric field 
~Edc is applied via the electrical gates and the STM tip. The 
resonance is driven by an additional small ac field. 
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ured to apply an electric field along the [001] axis. The 
STM tip serves as the third gate for spin manipulation, 
and as a contact for initialization and detection. Tak- 
ing advantage of the (110) natural cleavage plane (which 
lacks surface states), the applied electric field is confined 
in the (110) plane and the orientation is specified by the 
angle 9 from the [001] axis. 

An isolated Mn atom has a half-filled 3d shell and the 
spins of all five 3d electrons are aligned (Hund's rule) 
to form a S = 5/2 ground state. In GaAs a hole in 
the valence band compensates for the differing valences 
of Mn and Ga. We describe the core spin-valence hole 
dynamics with the following effective spin Hamiltonian: 

J^ pin = aS-s + /31-s, (1) 

where 1 and s are the orbital angular momentum (1 = 1) 
and the spin of the bound hole respectively. Our tight- 
binding calculations [24| estimate the exchange coupling 
a and the spin-orbit coupling j3 to be about 300 meV and 
—80 meV respectively. The exchange interaction binds 
the valence hole with spin antiparallel to the Mn core 
spin with a binding energy of 113 meV [2^]. The spin- 
orbit interaction in GaAs configures the orbital angular 
momentum of the hole parallel to its spin. The total 
angular momentum of the (Mn core + hole) complex is 
J = S+l+s, and the ground state of this complex has J = 
1 |2^| (both 1 and s are antiparallel to S), confirmed via 
ESR [26(. Our proposals for spin control involve energy 
scales smaller than a or (3, so only the lowest energy 
multiplet with J = 1 is of interest here. 

The degeneracy of the J = 1 Mn ion can be substan- 
tially split by external electric fields, and the eigenstates 
depend strongly on the electric field direction. This will 
be the source both of state splitting (analogous to the 
static magnetic field in traditional ESR) and state cou- 
pling (analogous to the oscillating perpendicular mag- 
netic field in traditional ESR). We find the following 
electric-field-dependent Hamiltonian: 

J^(E) = j[E x (J y J z + J z J y )+c.p.} , (2) 

where E is an electric field, c.p. stands for cyclic per- 
mutation, and {x, y, z} stand for the 3 major axes of the 
cubic crystal. Note that this Hamiltonian does not break 
time-reversal symmetry, for the angular momentum op- 
erators J always appear in pairs. We calculate, using the 
probability densities of the hole state found in our tight- 
binding calculations and first-order perturbation theory, 
7 = 6.4 x 10~ 30 Cm, corresponding to jE = 160 fj,eV 
for E = 40 kV/cm. This exceptionally large splitting is 
equivalent to that generated by applying a 1 Tesla mag- 
netic field using the measured g- factor (2(|, 2.77. The 
linear dependence on electric field, critical to producing a 
large splitting, originates from the lack of inversion sym- 
metry of the substituted ion in a tetrahedral host. The 
energy splittings from an electric field applied to bound 



states at inversion-symmetric sites in crystals, or elec- 
trons bound in atoms or ions in vacuum, would depend 
quadratically on the electric field and would be corre- 
spondingly much smaller. The other essential element 
causing this large splitting is the larg e (~ 10A) Bohr ra- 
dius of the bound valence hole [24ll27| Recent progress in 
theory and scanning tunneling microscopy of Mn dopants 
in III-V semiconductors has confirmed the large sp atial 
extent of the bound hole wavefunction |U 0, |2s| . 
Thus the response of the Mn wavefunction to electric 
fields is substantial compared to other ion levels associ- 
ated with transition- metal (magnetic) dopants. 

In the basis \X), \Y) and \Z), defined by J a \a) = 0, 
the Hamiltonian can be written as 

/ E z E y \ 
JT/(E) = -jE \ E z E x . (3) 
\E y E x / 

The energy eigenvalues in units of ^E are the roots of 
the characteristic polynomial, 

x 3 ~ x + 277 = , (4) 

where 77 = E x E y E z . A static electric field E^c splits all 
three eigenstates in energy except when the field is in 
the [111] direction (or equivalent), for which two of the 
eigenstates remain degenerate. 

The energies of the three states are £1 = (— cos 9 — 
\/4- 3 cos 2 6>)/2, £ 2 = (-cos 9 + V^-Scos 2 9)/2, and 
£3 = cos6>, shown by the solid, dashed, and dotted 
curves respectively in Fig. 0V). The eigenstate (£3) = 
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FIG. 2: (color online) Fhe ionic spin system as a function of 
the dc field orientation, (a) The energies of the J = 1 states £1 
(solid), £2 (dashed), and £3 (dotted), (b) The corresponding 
eigenvectors parametrized by the angle O. (c) The coupling 
between and ^2) due to the ac field, (d) The scaled LDOS 
of the two possible final states (solid) and \&) (dashed) 
probed four monoatomic layers directly above the Mn dopant. 
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(\X) — \Y))/i/2 is independent of 6. The independence 
of | £3) from E (in this geometry) motivates us to de- 
fine a pseudospin 1/2 constructed from the other two 
states, and l^)- These eigenstates can be writ- 
ten as |£i) = (sinO/V2, sinG/\/2, cosG) and \&) — 
( — cos G I v2, — cos G / v2 j sin G) , where G is the angle of 
between and the \Z) basis (Fig.^b)). Note that all 
the eigenvectors are real because of time-reversal symme- 
try. 

Preparation of the initial pseudospin state is achieved 
by applying an electric field to split the state energies, 
and allowing the hole to relax into the ground state. The 
electric field from the STM tip locally bends the bands of 
the semiconductor and permits ionization of the bound 
hole; this has been demonstrated for Mn in GaAs |27ll29| . 
Rapid initialization of a high purity pseudospin state can 
be achieved by using the local band bending effect to 
move the two higher-energy levels (£2, £3) to the position 
shown in Fig.|3{a), so a hole in those states would ionize 
and be replaced by a hole in the lowest energy state (|£i)). 
At a temperature of 0.5 K and a dc field of 100 kV/cm, 
the occupation of the next highest state (|£a)) would be 
less than 10 -4 . We have chosen Ed c such that 16*1 < 
(tt — tan -1 \/2), so that (not ^3)) is the ground state 
(see Fig. Eta)). Band bending also changes the effective 
radius of the bound hole wave function; gate voltages 
applied at the surface could thus control the coupling of 
two bound hole states in an analogous way to approaches 
111 Refs. 0[li for quantum dots and donor states. 

In order to manipulate the initialized spins the tip- 
sample bias should be increased adiabatically (slower 
than h/^Edc)) to bring all three levels into the semi- 
conductor's energy gap (see Fig. CJb)). This shift with 
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FIG. 3: (color online) Schematics of controlling the spin states 
via local band bending. The dot-dashed lines show the chem- 
ical potential. Shaded regions are filled states. CB and VB 
label the conduction and valence bands of the semiconductor, 
(a) Initialization. For this voltage occupation of the |£i) state 
dominates, (b) Manipulation: Bring all the states into the 
gap, but control the bias voltage below the threshold where 
the current starts to tunnel through these states. The oscil- 
lating field (E ac ) drives transitions between the and the 
^2} states, (c) Detection: Bring the final state further into 
the gap, so that electrons can tunnel from the tip into the 
acceptor state. The final state is identified according to the 
amplitude of the tunneling current (Fig. |2Jd) ) . 



bias is described for Mn in p-doped GaAs in Ref. [27} • 
The bias voltage has to be maintained below the critical 
value at which electrons start to tunnel directly through 
these levels, so that the transitions between these states 
remain coherent. Spin resonance can now be driven by 
applying a small oscillating electric field E ac (i) to the 
static field ~Edc- The Hamiltonian 

^esr = ^/(Edc) + J#i\Ei ac (t)] . (5) 

To have a well-defined pseudospin 1/2, constructed out of 
|£i) and \^), the coupling of these two states to (£3) must 
vanish. For the schematic in Fig. ^ the oscillating field 
can be applied either along the [110] direction through 
the STM tip or along the [001] direction through the 
gates. Both choices leave ^3) unaffected and only couple 
|£i) and (£2) to each other. To see how the states are 
coupled by the ac field, we write out J#j[E ac (t)] using 
the eigenstates of J^i(Edc) as bases. We assume that the 
ac field E oc (t) is along the [110] direction. 

/ -sin2G cos2G \ 
Jfi[Eac(t)] = jE ac co$(ujt) cos2G sin2G (6) 

V 0/ 

The off-diagonal term cos2G, plotted in Fig.^c), shows 
how the coupling between the two coupled states changes 
with the field orientation. The coupling is maximized 
when the static field is completely along the [001] direc- 
tion (8 — 0). Then the static and oscillating electric fields 
are perpendicular to each other, just as the static and os- 
cillating magnetic fields are perpendicular to each other 
in traditional ESR. In the limit E ac -C Edc, the diago- 
nal term can be neglected and our configuration works 
just like conventional ESR. The Rabi frequency obtained 
from the standard Rabi formula is 



nn = -y(jE ac cos2G) 2 + (hu - 7^ C V4 - cos 2 o) 2 . 

(7) 

For E ac = E dc /4 = 25 kV/cm, and G = tt/2, Q/2-k = 
12 GHz, corresponding to a Rabi time of 80 ps. Ensemble 
spin coherence times T 2 * measured by traditional ESR in 
GaMnAs exceed 0.5 ns (several times the estimated Rabi 
time), and appear due to the inhomogeneous environ- 
ments of Mn ions [2(J; the T^'s of individual spins are 
expected to be considerably longer. Hyperfine interac- 
tions, which significantly affect conduction electron spin 
coherence, are expected to be weak for Mn ions as the 
overlap of the valence p orbitals with the nucleus is small. 

High-fidelity determination of the orientation of the 
pseudospin can be achieved by measuring the total tun- 
neling current through the final state with the STM |IJ. 
When the tip-sample voltage is increased, and the semi- 
conductor bands bend further (see Fig. EJc)), current 
starts to tunnel through the bound hole wavefunction 
state |27l |29| and the tunneling current is proportional 
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to the probability density of the state at the STM tip 
location. The spatial structure of these J = 1 states 
is highly anisotropic H (23, H. We calculate the two 
eigenstates to have the following spatial structure, 

(r|&> = c i XY (r\X+Y)+c i z (r\Z) , (8) 

where c^y — sin 6, c? XY = —cos 9, and \X + Y) = 
(\X) + \Y))/V2. If the STM tip is positioned directly over 
the Mn dopant, it probes the nodal plane, (r\Z) = 0. In 
this particular case, the local density of states (LDOS) is 
simply proportional to the projection on to the (r\X + Y) 
state, which is |sin6| 2 for and |cos0| 2 for |£ 2 )- 

Thus for a static electric field along 9 = (O = tt/2) 
and the pseudospin in the |^) state no current will be 
detected, but for the pseudospin in the state there 
will be current detected. The difference in current for the 
|£i) and |^) states is shown in Fig. Etd). For this posi- 
tion about 10% of the LDOS is not spin dependent, which 
reduces maximum visibility to 90% (at, e.g., 6 » 0). Spa- 
tial averaging of the LDOS over a typical experimental 
2 A changes the visibility by only a few percent. The 
asymmetric angular dependence is due to the lack of in- 
version symmetry of the substituted ion in a tetrahedral 
host. Current measurement timescales can be very fast 
as STM exp eriments performed at 50 GHz have demon- 
strated [2fJ. We also assume that the tunneling current 
is small so that spin flip does not occur during the mea- 
surement. 

Controllable coupling of two spins permits use of these 
Mn ions for quantum information processing. Estimates 
of the overlap of holes bound to two separated Mn 
ions [24| indicate ~ 100 meV splittings of Mn pair states 
for ions separated by 12 A along the (110) direction. The 
overlap falls off for larger separations according to the 
- 13 A wave function radius of the bound hole 9], so 
it would be ~ 0.1 meV for two ions 10 nm apart. This 
overlap could be reduced, increased, or eliminated with 
a gate between the two ions mini. By using single-Mn 
manipulations to put single-ion quantum information in 
the proper pair of single-Mn states, the Mn pair state 
splitting can be used to perform CNOT operations in an 
analogous way to how the singlet-triplet splitting is used 
for a CNOT with spin- 1/2 qubits. 

In conclusion, we have presented a concrete proposal 
for electrically initializing, manipulating, and detecting 
single pseudospin states of a magnetic dopant in a semi- 
conductor. All-electrical spin manipulation should be 
possible for other impurities in tetrahedral semiconduc- 
tors characterized by J > 1/2 ground state spins (e.g. 
most transition metal ions in most tetrahedral semicon- 
ductors, or the nitrogen- vacancy center in diamond). In 
a future scalable architecture the STM tip would be re- 
placed by a gate-controlled contact. The controlled resis- 
tance of that contact would permit alternation between 
the gate configuration for manipulation and the contact 
configuration for initialization and detection, all without 



moving parts. The [001] static electric field, here assumed 
to be implemented with gates, may also be replaced by 
an internal electric field from a doping gradient (such as 
inap-n junction) , or even a static strain field. The Mn 
ions could be controllably placed within the surface rel- 
ative to the contacts using current pulses from an STM 
tip as described in Ref. |23j. 

This work is supported by ARO MURI DAAD19-01- 
1-0541 and DARPA QuIST DAAD-19-01-1-0650. 
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